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ABSTRACT 

A  similarity  is  found  betwee^^  the  behavior  of  certain  waves 
in  a  plasma  in  a  magnetic  field  and  gravity  waves  on  the  svirface 
of  water  of  finite  depthj  in  both  cases,  the  particular  motion 
studied  in  this  report  is  caused  by  an  uniformly  moving  piston, 
started  impulsively  from  resto   This  similarity  in  behavior 
develops  asymptotically  in  time  away  from  the  initial  state,  both 
in  the  linearized  description  of  the  motion  and  in  a  more  general 
non-linear  small-pert\irbation  description  which  includes  the 
linear  asymptotic  motion.   The  approximate  non-linear  theory 
yields  the  third  order  differential  equation  2w„  +  3wWj-  +  w^-j-j-  =  0 


which  can  be  interpreted  as  describing  a  reversible  dispersion 
process,  in  contrast  to  Burgers'  well-known  second  order 
differential  equation  which  describes  an  irreversible  diffusion 
pro cess o 
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NYO-9082 

SIMILARITY  IN  THE  ASYMPTOTIC  BEHAVIOR  OP  COLLISION-FREE 
HYDROM^iGNETIC  WAVES  AND  WATER  WAVES 

1„   INTRODUCTION 

In  the  main  body  of  this  report,  we  study  the  motion  of  a  collision- 
free  plasma  in  a  magnetic  field  caused  by  a  piston^  started  impulsively 
from  rest,  moving  with  a  small  uniform  velocityo   In  particular  we  are 
Interested  in  the  asymptotic  time=dependent  behavior  near  the  equilibrium 
stateo   Since  the  study  of  the  corresponding  motion  of  water  waves  is  so 
similar,  the  description  is  given  in  the  Appendices,,   An  interesting 
difference,  however,  is  that  the  original  water  wave  motion  is  two- 
dimensional  and  becomes  quasi-one-dimenslonal  asymptotically  for  large  time, 

We  consider  an  idealized  fluid  model"  of  singly-charged  particles 
moving  in  an  electromagnetic  field.   The  main  properties  of  this  model  are; 
i)   Charge  neutrality  ~-  the  basic  property  of  a  plasma, 
ii)   No  collisions  --  the  charged  particles  interact  only 
through  the  self-consistent  electromagnetic  fieldp  and 
iii)  Zero-temperature  plasma  =>-  completely-ordered  motion,  in 
which  the  fluid  motion  is  the  same  as  the  particle  motiono 
The  physical  implications  of  these  assiiraptions  and  a  study,  primarily  of 
equilibrixom  solutions,  are  described  by  Adlam  and  Allen""  ( 1958a )j  Davis, 
Lust  and  Schluter  (1958)  a'^d  others.""" 


Some  properties  and  solutions  of  this  model  have  been  summarized  by 
Gardner  et  al  at  the  Geneva  Conference  (1958)o 

'Also  some  related  calculations  describing  behavior  near  the  Initial 
state  of  a  pinched  discharge  are  given  in  a  contribution  to  the 
Geneva  Conference  ( 1958b  )« 


"Eogs,  Ao  Banos  and  M,  H„  Mittelman  at  UCRL,  Livermore,  and 
M„  Do  Kruskal  at  Princeton  (private  communication) o 


k 


The  study  of  water  waves  Is  an  old  discipline  aid  the  differential 
equations  governing  their  motion  have  been  known  for  many  years. 
Actually  however,  the  difficult  mathematical  analysis  of  basic  non- 
linear free-boundary  type  problem  is  still  far  from  a  fully-developed 
state,  eogo  see  the  treatise  by  Stoker  (1957)o   The  formulation  for 
water  of  finite  depth  is  s-uramarized  in  Appendix  II  where  the  non-linear 
approximation  is  carried  outc 

The  motion  which  we  study  of  the  plasma  and  water  waves  may  be 
only  roughly  similar  near  the  initial  stated  and  the  analogy  between 
the  two  phenomena  develops  only  after  some  time  has  passed  and  the 
motion  is  approaching  the  equilibrium  state.   This  similarity  in 
behavior  exists  asymptotically  in  the  linearized  description  (  8   3  and 
Appendix  I),  as  well  as  in  a  time=dependent  non-linear  approximation 
(  6  i\.   and  Appendix  II)  which  includes  both  the  linearized  asymptotic 
motion  and  the  final  equilibriijm  state.   The  range  of  validity  of 
linearized  theory  and  the  approximate  non-linear  theory  is  indicated 
graphically: 

-^  ____.^=_=_...^ Full  non-linear  problem 

i, Linear  theory 


-> 


-Approximate  non-linear  theory- 


Initial 
State 


Equilibrium 
State 


->  t 


2«   FORMULATION 

As  a  consequence  of  the  basic  property  of  charge  neutrality,  the 
motion  of  a  cold  plasma  consisting  of  singly-charged  particles  can  be 


described  by  a  one-fluid  theory.   The  state  of  the  plasma  is 

characterized  by  a  system  of  five  first-order  partial  differential 

equations  for  the  variables  (n,  u,  v,  B,  E)  where 

n  =  n,  =  n  ,  number  density  of  ions  (or  electrons) 

u  =  u  =  u  p  x-component  of  the  ion  (or  electron)  velocity 

V  =  (v  "V  )p  difference  between  the  y-component  of  ion  and 

electron  velocities 

B,  z-component  of  the  magnetic  field 

Ej,  y-component  of  the  electric  field 

ix) 
and,  with  the  exception  of  E  \    x-component  of  the  electric  field 

(charge  separation  field),  all  other  components  are  zeroo   These  six 
quantities  depend  on  the  single  space  coordinate  x  and  the  time  to   For 
convenience  we  consider  these  quantities  to  be  normalized  in  the  follow- 
ing way: 

m      n 
o 

u:   A  =  [B^/nn„(m,  +  m  )]'^^^ 
o     00+    = 

?  T  /? 

v:   a  =  [B  (m,  +  m  )/un  ra.m  ]  ' 

o      0+     -'■^0  +  =' 

B;   B 
o 

E:   E  =  A  B^ 
coo 


E 


(^^  e(^^=  a  B 


o  o 


xs  X  =  [m,m  /ue,  n^(m,+m  )] 

t:  t^  =  [m^mye/B/]^/2  =  x^A^ 

where  \i   is  the  magnetic  permeability  of  space,  e^  and  m^  are  the  charge 
and  mass  of  the  ion  (minus  sign  for  electrons),  and  the  subscript 
(zero)  refers  to  a  constant  initial  or  boundary  statco   Then  the  re- 
sulting normalized  equations  are  one  conservation  eauation  of  number 


density  and   two   of  momentxim  and  Maxwell's   equations    (two): 

n^   +    (nu)      =  0  (1) 

t  X 

u.    +   uu     =  vB  (2) 

t  X 

v^   +  uv     =  E   -   uB  ( 3 ) 

t  X 

B.    +   E     =   0  ik) 


B      =   -nv  (5) 

X 


where  the  displacement  current  is  neglected  in  (5)  and  lettered  sub- 
scripts indicate  partial  diff erentiationo  In  this  one-fluid  theory, 
the  x-component  of  the  electric  field  is 


/  \     iti,  =  III 

E^^^  =  (-i,  Tr^)vB  (6) 

m,  +  m^ 


Also  "by  eliminating  the  y-component  of  the  electric  field  E  between 
(3)  and  (Lj.)  we  get  a  simplificatlono   Then 


(v^  +  B)^  +  [u(v^  +  B)]^  =  0  (7) 


Comparing  (7)  with  (1),  (v  +  B)  satisfies  the  same  equation  as  n, 
Thus  we  can  set 


V  =  n  -  B  (8) 

X 
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If  this  is  true  Initiallyo   Now  we  have  a  reduced,  but  equivalentj, 
system  of  fovir   equations  (l),  (2),  (5)  and  (8)  for  (n,u,Vj,B)  in  place 
of  the  original  five  and  (3)  may  be  regarded  as  the  defining  equation 
for  Eo   We  note  that  the  particle  path  and  the  plane  t  =  constant  are 
the  only  set  of  real  characteristics  (actually,  a  double  set)  for  this 
differential  equation  systemj  and  neither  A  nor  a  are  characteristic 

'^      '  GO 

velocities  in  this  senseo 

For  later  use  {f   l|)  we  look  for  equilibrium  solutions  of  the 
reduced  system  (1),  (2),  (5)  and  (8),  noting  that  the  equations  are 
Gallilean  invariant,  loeo,  invariant  under  the  transformations 
jx   — >  (x-u  t) ,  t  — s»t,  u  — *  (u~u  )l  where  u  is  constanto   Setting 

a 

mass,  momentum  and  energy  follow  by  inspection 


—  =  0  in  (1)  and  (2),  combined  with  {$)  and  (8),  three  Integrals  of 


nu  =  u 
o 


(9) 


o     2     o    2 

u  2 
|(u^  +  v^)  +  B  =  -|-  +  B.  (11) 

Combiniag  (10)  and  (11)  yields  the  "phase  relation" 

v2  =  (B-l)2[l  -  (B±l)!  ]  (12) 

Eliminating  n,  u,  and  v  by  (5),  (9),  (10)  and  (12),  3(x)  can  be 
integrated  explicitly  ( cf «  references  in  ^  l)o   We  obtain  the  solution 
in  the  limiting  case  of  a  small  disturbance  by  the  transformation 
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4  =  e^/^x  and  letting  n  =  1  +  en^^^C),  u  =  1  +  eu^^^C),  v  =  e^/^v^-^^C), 

B  =  1  +  eB   (C)  and  u  =  1  +  ea,  where  a  =  const,  to  bo  determined  and 

e  =  (B   -1)«  1  is  the  magnetic  field  amplitude.   Then  (5),  (9),  (10). 
max 

and  (12)  become 


#i^=-v'l'  (13) 

dt, 

n(l)  +  u^^^  =  a                                     ilk) 

u(l)  +  b(1)  =  a                                     (15) 


(v(l))2=  (B(^^)2(2a  -B^^h  (16) 


Eliminating  v^-""^  between  (13)  and  (l6),  B^-^^  satisfies 


dB!£=  _3(l)(2^„3(l))l/2  (   ) 

dc 


which  yields  the  equllibri'um  solution 

B^^^  =  2a  sech^(a?)  (l8) 

We  choose  a  =  l/2  so  that  B    =1,  consistent  with  our  definition  of 
'  max    ' 

eo   (l8)  is  a  pulse^type  solution  similar  to  the  small- amplitude 
solitary  water  wave,  eogo  see  Keller  (19U.8),   The  boundary  conditions 
chosen  in  (9),  (10)  and  (11)  have  eliminated  some  additional  equilibrivmi 
solutions  (periodic)  which  can  be  recovered  from  the  time-dependent 
non-linear  approximation  derived  in  ^  I).© 
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3,   LINEAR  APPROXIMATION* 

To  describe  the  motion  caused  by  small  initial  distiorbances  we 
linearize  the  reduced  plasma  equations  (1),  (2),  (5)  and  (8)  by 
perturbing  on  the  rest  state  (npU,v,B)  =  (1,0,0,1) o   The  resulting 
equations  in  the  first  approximation  are 


°t'^'  +  u^'^'  =  °  <19) 


u^'l'  =  V<1'  (20) 

B  (1)  =  =v'l'  (21) 


V  (1)  =  n<l).  B<1»  (22) 


By  elimination  we  obtain  a  single  fourth  order  differential  equation 

(1) 
for  the  x-component  of  the  velocity  u    t 


Also  n   p  B'    and  v    satisfy  the  same  equation  (23)o   We  now 
consider  the  plasma  motion  for  x  >  0  caused  by  an  electrically- 
neutral  piston""  (nominally  at  x  =  0)  started  impulsively  from  rest 


The  results  in  this  section  were  presented  by  Gardier  (1958)  at  the 
Controlled  Thermonuclear  Conference  in  Wash,  Dc  Co  and  have  been 
summarized  in  the  Proceedings  of  the  Geneva  Conference  by  Gardier  et  al 
(1958)o 

.  This  idealized  problem  is  considered  for  convenieice.   We  could  just 
as  easily  take  the  same  differeitial  equation  (23)  satisfied  by  B^-'-^ 
and  consider  a  magnetically=driven  pistono 
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and  moving  with  uniform  velocity  U^  in  the  positive  x  direction.   Then 
the  initial  conditions  are 

u^^^x,0)  =  u^^^^x,0)  =  0  (2I4.) 

and  the  boundary  conditions  for  t  >  0  is 

^^^\(Q>^)  =  H(t)  (25) 

1 

where  H(t)  is  the  unit  step  function  (=  0  for  x  <  0  and  =  1  for 
X  >  0)0   Equation  (23)  with  the  conditions  (2l|.)  and  (25)  can  be  solved 
by  applying  the  Laplace  transform  with  respect  to  the  time  t.   The 
solution  is 

n^lu^^^      r  ^,^p  [ts[l-a(l+s2)-l/2]]      (26) 
1  I 

where  a  =  x  /t  and  the  contour  /^  goes  from-iooto  +  iooin  the  right 
half  of  the  complex  s-plane.   To  carry  out  the  inverse  Laplace 
transform  (26)  appears  difficult"  but  the  asymptotic  behavior  for  large 
t  can  be  determined  by  the  method  of  steepest  descent.   Then  the 
asymptotic  form  of  the  velocity  distribution  for  large  t  (and  x)  is 


u^^^x.t)  ^ 


00 


Ai(«^)do6  (27) 


where  p  =  (2/3)-^^^(x-t )/t^/^  and 


-  11 


Ai(p)  =^ 


oo        .3 

dX   cos(-^  +  BX)  (28) 


is  the  Airy  function.   The  same  asymptotic  behavior  is  found  'for  the 
linearized  motion  of  water  waves  in  Appendix  I,  equation  l-17o   Equation 
(27)  describes  a  shock- type  velocity  distribution  with  damped  oscilla- 
tions behind  the  "shock"  front  ((3  <  0)o   The  velocity  decreases 
exponentially  ahead  of  the  front  P  >  0,  and  for  negative  p  the  velocity 

oscillates  about  and  approaches  the  piston  velocity  U^  as  p  — *■     -00 » 

-3/1 
Also  the  amplitude  of  the  oscillation  decreases  like  p"  '  ^  as  P  -*=•  -cr  « 

However  this  wave  is  not  a  shock  in  the  usual  sense  since  it  is  non- 

1/3 

stationary,  widening  with  time  like  t  ^ -^l    thus  (2?)  describes  a  slowly 

dispersing  waveo 

l|o   TO^r-LINEAR  A'^PROXIMATION 

So  far  the  only  apparent  relationship  between  the  linearized 
asymptotic  behavior  of  plasma  motion  (equation  (2?)  in  ^  3)  and  the 

pulse-like  equilibrium  solution  (equation  (I8)  in  o  2)  is  their 

?  1/2 

propagation  speed  A  =  [B  /\in    (m^  +  m_)]  ^      as  e  — >  Oc   In  this 

section  we  seek  that  approximate  time-dependent  description  which  is 

valid  for  large  time,  ioeo,  valid  under  the  conditions  that  it 

1)  includes  the  equilibrium  solution  (I8)  and  2)  is  in  the  same 

asymptotic  direction  in  the  x-t  plane  as  the  linearized  asyinptotic 

1/3 
solution  (27),  that  is,  the  direction  defined  by  P  =  (x-t)/t  ^   = 

const o   It  is  clear  from  (17)  and  (I8)  that  a  non-linear  approxima- 
tion is  required.   A  comparable  formulation  has  been  fouid  in  the 
approximate  description  of  the  motion  of  flood  waves  in  rivers  by 
Morikawa  (1957). 
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To  obtain  the  desired  non-linear  approximation  of  the  plasma 
equations  (1),  (2),  (5)  and  (8),  we  first  apply  the  above  condition 
1)  which  states  that  the  pulse  solution  (l8)  is  to  be  included  in 
our  approximate  formulation.   Thus,  we  make  the  following  scale 
transformations: 


?  =  e^/^(x-t)  (29) 


and 


V  =  e^/^Vo  (30) 


Condition  2)  above  fixes  the  necessary  scale  transformation  of  time: 

l/3 
since  we  want  (x-t)/t  '       to  be  invariant  with  respect  to  the 

coordinate  transformations  (29)  and  (31).   In  addition,  we  make  a 

perturbation  expansion  with  respect  to  e  on  the  rest  state  (n,u,v,B)  = 

(1,0,0,1): 


n  =  1  +  en^^^?,r)  +  t^n^^\^,t)    +  «,.  (32) 

u  =  eu^^^?,r)  +  e2u^2^5,r)  +  oc.  (33) 


V  =  eV^^^U,r)  +  e^V^2^?,r)  +  ...  (31;) 


B  =  1  +  eB^^^cr)  +  e^B^2^^,r)  +  ...  (35) 
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By  the  scale  transformations  (29),  (30)  and  (31),  the  equations 
(1),  (2),  (5)  aad  (8)  become 


er^  +  [n(u-l)]^  =  0  (36) 

ev^   +  (u-l)u^  =  0  VS  (37) 

B^  =  -nV  (38) 

eV^  =  n-B  (39) 


Putting  (32  to  35)  in  (36  to  39)  and  equating  the  coefficients  of 
powers  of  e  in  the  usual  way,  we  obtain  the  equations  which  the  lowest 
order  approximation  (n^-*"^,  u^"""^,  V^-^%  B^""-^)  satisfies: 


n(l)  -  u^^)  =  0  (1+0) 

u^l)  =  b(1^  =  -V^l)  (1,1) 

n^l)  ^   B^^^  =  0  (1,2) 


and 


„<2>  -  b(2)  =  v(l)  (U6) 

(2)     (2)     (2)     (2) 

The  second  order  approximation  terms  (n    ,u    ,V    ,B    )  may  be 
eliminated  among  (1,3  to  I;6)  and  u^"""^  (or  n^"*"^  or  B^^' )    satisfies  the 
third  order  non-linear  equation 
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2u 


(1) 


■r^ 


+  3u^^^4^^  +  4U  =0  (1^7) 


(14.7)  describes  the  plasma  motion  traveling  at  the  velocity 

A  =  [B^  /p-n^dn,  +  m  )]  '   in  the  positive  x  direction  near  aid  in- 
o     00+    — 

eluding  the  equilibriiom  state.   An  equation  of  the  same  form  describes 
the  motion  of  water  waves  traveling  at  the  speed  of  small  gravity 
waves  (Appendix  II,  equation  (11-32) )» 

In  contrast  to  Burgers'  second  order  equation  which  describes  an 
irreversible  diffusion  process,  ([<.7)  describes  a  dispersior  process 
which  is  reversible  in  the  following  sense:   the  equation  is  invariant 
with  respect  to  change  of  sign  £,   — *•  -£,   and  f — >  -to      The  solution 
of  the  initial  value  problem  with  a  step-function  initial  distribution 
u^-^hz,0)    =  I  -  H(^)  (H  =  0  for  C  <  0  and  H  =  1  for  ?  >  0)  would  be  of 
some  interest.   But,  the  Hopf  (1950)-  Cole  (l95l)  technique,  which 
yields  the  corresponding  solution  of  Biorgers*  equation,  apparently  does 
not  apply  here.   We  plan  to  carry  out  a  numerical  Integrationo 

It  is  easily  shown  how  the  linearized  asymptotic  solution  (27)  is 
related  to  (l|-7)  and  that  the  pulse  solutio-^  (I8)  is  included  in  the 
equilibrium  solutiois  of  (14-7).   Neglect!  ig  the  non-linear  terra  in 
( I4.7 )  we  look  for  similarity  solutions  of 

2uj.^^  +  u^l^  =  0  (1|8) 

of  the  form  u^-""^  ='C"^f(8)  where  p  =  b^t^  and  b,  m  and  n  are  constants o 
Then,  m  is  arbitrary  and  n  =  -l/3  and  f  satisfies  the  ordinary 
differential  equation 


Burgers'  equation  is  w.+  w  w^  =  v  w   where  v  =  const. 
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f  "  '  -  pf '  +  3mf  =  0  ( 14.9 ) 


1/3 
where  b  =  (2/3)    for  convenlenceo   We  consider  three  particular 

values  of  m  =  0,  -l/3,  -2/3: 

Case  1,  m  =  Oc   For  this  case  (  1_l9  )  becomes 


(f')"  -  p(f')  =  0  (50) 


and  f   satisfies  the  Airy  equation  so  that 


f  =  A^(P)  (51) 

The  linearized  asymptotic  solution  (2?)  corresponds  to  this  case^ 
Case  2,  m  =  -l/3«   For  this  case  a  solution  of  ( li? )  is 


f  =  A^(P)  (52) 

and 

Jl)  ^^-1/3  ^^(0)  (53) 


(53)  is  the  asymptotic  solution  obtained  by  Jeffreys  aid  Jeffreys 
(1956)  in  their  study  of  the  motion  of  water  waves. 
Case  3,  m  =  -2/3o   For  this  case  a  solution  of  (14-9)  is 


f  =  A^'(S)  (514-) 
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(53)  is  unbounded  for  p  — >  -oobut  this  case  is  of  interest  since 
the  non-linear  equation  (I4.7)  also  has  a  similarity  solution  of  the 
form  u^"''^  =T"^/-^f(p)  where  f(p)  satisfies  the  equation 


111     I     t  /  ^^» 

f   +  ff  -  pf  -2f  =  0  (55) 


It  would  be  of  interest  to  determine  whether  there  is  a  bounded 
solution  of  (55)  over  the  entire  range  of  -co    5  P  £  00.  Equilibrium 
solutions  of  (I4.7)  are  obtained  by  seeking  solutions  of  the  form 
u  '    =   F(4-aT)  where  a  =  const o  to  be  determined.   Since 
(?-af)  =  e^[x   -    (l+ea)t],  (ea)  is  the  velocity  deviation  of 
the  steady  progressing  wave  from  the  critical  speed  A  «   By  ( I4.7 ) , 
P(^  -  slT)    satisfies  the  equation 


P'"  +  3Pf'  -  2aP  =  0  (56) 


I 

Multiplying  (56)  by  P  and  integrating  twice  we  obtain 


(p')2  +  p3  _  2ap2  =  c^P  +  c^  (5?) 


where  c,  and  Cp  are  integration  constants.   If  we  Impose  the  condition 

t    1 1 

that  the  solution  vanishes  at  Infinity,  then  F  =  p  =  p  =0  Implies 

that  c,  =  Cp  =  0°  and  (5?)  becomes 

(p')2  =  p2(2a  -  p)  (58) 
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which  corresponds  to  the  solitary  wave  equation  ( 1? ) »   If  we  define 
e  as  the  wave  amplitude,  then  uA^  =  1  implies  that  a  =  l/2o   The 
flow  velocity"  is  supercritical  at  infinity  and  subcritical  at  the 
peako   Periodic  solutions  are  obtained  by  relaxing  the  conditions 
at  infinity o   For  example  P  =  P  =0  implies  that  Cp  =  0  and  we  obtain 
for  fixed  amplitude  a  one=parameter  (wave  length)  family  of  periodic 
solutions  ( cnoidal  waves)  with  speeds  less  than  or  equal  to  the 
solitary  wave  speedo 

Finally,  we  summarize  the  correspondeice  between  variables 
describing  the  asymptotic  motion  of  plasma  waves  and  water  waveso 

Plasma  Water 

A„=  [  B//^n„  { m^4™_ )  ]  VS  ( gh  J  l/2 

3to=lm+m./M.e+ng(m^.+m_)l  h^ 

We  point  out  that  for  plasma  motion  the  y-component  of  velocity 

(1) 
V    is  normalized  with  respect  to  a  different  velocity 

P  T  /p         ( "]  ) 

a  =  [B   (m.+  m  )/an  m^m  ]  '   than  u    which  is  normalized  with 
o     o   +   =0  +  =° 

respect  to  A  (cfo  $   2) „      Thus  the  actual  y-component  of  velocity 

1/2 
exceeds  the  x-component  of  velocity  by  a  factor  a  /a  =  (my^m^)  '  o 

[1  +(m  /m.  )]  for  (mVm  )  »  1, 


Note  that  the  flow  velocity  is  from  left  to  right  with  respect  to 
the  solitary  wave  In  S  2   and  in  the  opposite  direction  in  this 
section,  appropriate  for  a  piston  moving  from  left  to  righto 
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APPENDIX  lo   LINEAR  APPROXIMATION  OP  MOTION  OP  WATER  WAVES 

To  compare  with  the  plasma  motion  described  in  6   2,  we  consider 
the  motion  caused  by  a  piston  (started  impulsively)  moviig  at  a  small 
uniform  speed  U^  into  water  initially  at  rest  of  \anlform  depth  h  » 
The  piston  (nominally  at  x  =0)  moves  in  the  positive  x  direction 
and  the  gravitational  acceleration  g  is  in  the  negative  y  directiono 
Neglecting  viscosity,  the  motion  is  described  by  the  velocity  potential 
0(x,y,t)o   Por  convenience  we  introduce  the  following  normalization: 

xs   h 
o 

tt        (hyg)l/2 

0t        (gh^)^/2 
\'      (gh^)^/2 

Since  the  boundary  conditions  on  the  water  surface  (cfo  Appendix  II) 
are  linearlzedj,  the  flow  is  nominally  confined  to  move  in  a  semi- 
infinite  strip  X  >  Oj  with  the  horizontal  bottom  at  y  =  0  and  the 
free  surface  at  y  =  lo   The  motion  is  goveried  by  the  Laplace  equation 
in  the  interiors 

^xx-^^y==0  (I-l) 

The  initial  conditions  are 

i2f(x,y,0)  =  0^(x,y,O)  =  Oo  (1-2) 

The  boundary  conditions  are;   at  x  =  0 


0^(0, y,t) 


X 


^1 


=H(t)  (1-3) 
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where  H(t)  is  the  unit  step  fiinctionj  on  the  bottom  y  =  0 

and  at  the  linearized  free  siorface  7  =  1 

The  solution  of  (l-l),  with  the  conditions  (1-2  to  5),  is  readily 
obtained  by  applying  the  Laplace  transform  with  respect  to  t  and  the 
Fourier  cosine  transform  with  respect  to  xj  the  resulting  second 
order  ordinary  differential  equation  in  the  remaining  variable  y  is 
easily  solvedo   After  a  simple  inversion  of  the  Laplace  transform, 
the  inverse  Fourier  transform  yields  the  integral  solution 

00 

0(x,y.t)=  =  ^         ^     cosCKa)   h^   ^gosh^^    cos[/<tanh  <)^/^tli  (1-6) 

The  wave  profile  y|(x,t)  =  0^(x,l,t)  and  from  (1-6) 

00 

QTT     it 

7)(x,t)  =  -   -^  ]    -'y|(Ktanh  )^)^/2cos(j(x)sin[)<tanhi()-'-/^t]      (1-7) 
We  make  an  asymptotic  evaluation  of  (1-7)  for  large  t  (and  x)  by  a 


variant  of  the  method  of  steepest  descent^   Since  the  integrand  of 

OD 

o 


f    1 
(I-7 )  is  regular  at)(^=  0  and  is  an  even  function  of  X»  then  1   =  ^ 

We  write  (1-7)  as  the  sum  of  two  sine  integrals 


00 
CD    ^  -^  -co 


cos  a  sin  b  =  -  i[sln(a=b)  -  sin(a+b)]         (I-8) 
1/2. 


where  a  =  X^,  b  =  (K.tanh>()^  t  and 


20  - 


(a+b)    =Xt[a  +    (tanhXy-^/^/X.-^/^]  (1-9) 

where   a  =  x/to      For  )(«   1    (large   x),    (1-9)    is   approximately 

.2 
(a+b)»Xt[a  +    (1-  g--)]  (I-IO) 

For  large  t  and  a  =  0(1)  (x>0)j,  the  integral  containing  sin(a-b) 

l/3 
dominates  that  with  sin(a+b)o   By  transforming  "}^=  (2/t)  '    \    and 

(a-1)  =  p/(2t^)-^/^j,  we  rewrite  (a-b)  in  (I-IO) 

\^        - 
(a-b)^7^fA-  +  e) 

and  the  asymptotic  solution  of  (1-7)  for  large  x  and  t  is 

U   r°°         3 
7|_(x,t)^2i     ^  sin(^  +  PX)  (l-ll) 

^oo 

Since  the  derivative  with  respect  to  p  of  the  integral  (l-ll)  is 

the  Airy  function 

oo 

Ai(p)  =  ^  I    dX  cos(^  +  pX)  (1-12) 

we  rewrite  (I-ll)  in  the  concise  form 


AGO 

(x,t)  ^U^J    Ai(Y)dY  (1-13) 

F 

where 

6  =  2^/^(x-t)/t^/^  (I=li|) 


• 


21 


The  wave  profile  (1-13)  is  of  the  same  form  as  (2?)  in  6  2o 
Comparing  (I-li4-)  with  (28),  p  differs  from  p  by  a  factor  3"^/^ 
(cf  Appendix  II,  equation  (II-32 ))  o   A  more  direct  comparison  with 
(27)  is  found  by  calculating  the  x-component  of  velocity 
u(x,l,t)  =  0  (x,l,t)  at  the  free  STorface  y  =  lo   From  (1-6) 


^^g^^^^^  =  ij       ^  sin(Xx)  I  l-cos[(x:tanhX)^/^t]j  (l-l5) 


Since 


J     "^ 

J  -00 

(I-l5)  becomes 


-^^  sin(K'c)  =  Ti,  X  >  0, 


^00 


^^g^^^^^  =1-1      ^     slnCXx)    cos[XtanhX)^^^t] 


(1-16) 


which  differs  from  the  wave  profile  (1-7)'  but  the  asymptotic  form 
of  (1-16)  is  the  same  as  (1-13)  and 


u(x,l,t)  ^^      \ 

"i  J 


00 


Ai(Y)dY  (1-17) 
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APPENDIX  II o   ASYMPTOTIC  NON-LI ^TEAR  APPROXIMATION  FOR  MOTION  OF  WATER 
WAVES 

In  two  diraensions  the  system  of  equations  governing  incompressible 
flow  are  the  conservation  equations  of  mass  and  mementim: 

u  +  V  =  0  (Il-l) 

X    7 


g  +  p  +1=0  (II-3) 

dt   ^y 

where 

The  variables  have  been  normalized  as  follows: 

t.  (hyg)^/2 

ps  pgh^ 

where  h  is  the  depth  for  the  rest  state,  g  is  the  gravitational 
acceleration  and  p  is  the  density  ( constant )«   In  addition  to 
(II-ls,2j,3)  we  impose  the  irrotatlonality  condition 

v  -  u  =  0  (Il-l;) 

X    y 

The  boundary  condition  at  the  horizontal  bottom  is 

v(x,o,t)  =  o;  (II-5) 
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and  at  the  free  surface 

v(x,h,t)  =  ^+u(x,h,t)^l  h(x,t)  (II-6) 

and 

p(x5h,t)  =  Oo  (II-7) 

Just  as  for  the  plasma  equations  (l),  (2),  (5)  and  (8),  the  particle 
path  and  the  plane  t  =  constant  are  each  double  characteristics  of 
(Il-l,  2,    3)o 

Following  the  same  procedure  by  which  we  derived  the  non-linear 
approximation  of  the  plasma  equations  in  c  L|.,  we  obtain  the  asymptotic 
approximation  for  water  waves  by  1)  including  the  small  amplitude 
solitary  wave  as  an  equilibrium  solution  and  2)  scali^g  the  time  so 
that  (x-t)/t  '   is  invariant,  as  dictated  by  the  linearized  asymptotic 
solution  (l-17)o   Then  the  required  scale  transformations  are  the  same 
as  those  for  the  plasma  equations  (cf<,  (29),  (30)  and  (31)) 

4  =  e^/2(x-t)  (II-8) 

r=  s3/2^  (,,.g) 

V  =  e^/^v  (11-10) 

and  the  remaining  variables  y,u  and  p  are  unchanged »   The  flow  equations 
(II-l,2j,3)  are  transformed  by  (11-8,9,10)  into 

u^  +  V  =0  (Il-ll) 

§  ^  P^  =  0  (11-12) 

e  11  +  py  +  1  =  0  (11-13) 

-  21;- 


where  the  notation  d(  )/d.'t  means 

^^  =  e(  )^  +  (u-1)  (  )^  +  Vo  (  )y  (II-IL^) 

and  the  irrotationality  condition  (ll-l;)  becomes 

eV_  "  u  =  Oo  (II-15) 

€    y 

In  addition,  the  boundary  conditions  (11-5, 6,7)  become 

V(C,0,r)  =  0  (11-16) 

V(?,h,t)  =  e  ^  +  [ti(€,h,r)-  1]  ^   h(^,T)  (11-17) 

and 

p(4,h,r)  =  0  (II=18) 

There  is  a  close  relationship  between  the  scale  transformations 
(ll"8,9,10)  and  the  transformation  of  variables  yielding  the  well-known 
long  wave  (or  shallow  water)  approximation  [cfo  Priedrichs  (19I4-8)  or 


Morikawa  (1956)] o   Both  transformations  yield  the  hydrostatic  approxi- 

dV 

y 

Morikawa' s  version  (differs  from  Priedrichs"  in  the  scaling  of  the 


mation  for  e  -«►  Oj,  ioe<,  in  (11-13)  ^  "a^  ^^    (p„  "*■  1)°   Also  in 


vertical  velocity  component)  the  Irrotionality  condition  (II-15)  states 
that  u  =  0  for  e  =  0,  ioeo  the  horizontal  velocity  component  u(4yT) 
is  independent  of  the  vertical  coordinate  yo   The  scale  transformation 


yielding  the  long  wave  approximation  differs  only  in  the  time  scaling, 

1/2 

being T=  e  '  t  Instead  of  (II-9)o   That  is,  the  non-linear  approxlma- 

3/2 

tionj  which  we  describe  in  this  Appendix  iX=   e    t),  describes  the 
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motion  at  later  time  and  closer  to  the  final  equilibrium  state  than 
the  long  wave  approximationo 

In  addition  to  the  scale  transformation  we  make  a  perturbation 
expansion  (cf<,  6  1|)  on  the  rest  state  (h,u,V,p)  =  ( 1,0,0,  (1-y)), 
Then, 

h  =  1  +  eh^^^?,r)  +  e^h^^^C,r)  +  ,oo  (11-19) 

u=  eu^^^?,y,r)  +  e2u^^^?,7,<t)  +  o.o  (11-20) 

V=  eV^^U^,7,t)    +  t^'J^^U^.j.r)    +    .00  (11-21) 

P  =  (l^y)  +  e[p^^^^,y,'r)  +  h^^^C,t)]  +  e^ip^^h^^jfi^)   + 

+   h^^^C,r)]  +  coo         (11-22) 

noting  that  at  the  free  surface  y  =  h 

v(?,h,r)  =  ev^i^?,i,r)  +  e^iY^^U^A.^)   +  v|^^4,i,r).h^^^C,r)]+ooo 

(11-23) 

p{^,h,t)    =   ep^^^C,l,T)+  e^[p^^h^>l»^)+   p|^^C,l,T)°h^^^?,T)]+  ooo 

(II-2l|) 

Putting  (11-19  to  22)  in  (ll-ll  to  l8)  and  equating  coefficients  of 
like  powers  of  e,  the  lowest  and,  if  desired,  higher  order  approxima- 
tions can  be  calculated  in  a  systematic  wayo   We  summarize  the 
equations  yielding  the  lowest  order  approximation  (h   ,  u   ,  V    , 
p^-^^)o   Then,  u^-'-^(^,T)  is  independent  of  y,  p^-'"^eO  and 

h^l)  »  J^^    =   0  (11-25) 

v^l^  =  -y  u^^^  o  (11-26) 
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AlsOj, 


u^^^  =  -(7^2)41^  +  P(4,r)  (11-27) 

V^2)  =  (y3/6)u^l^  .  yp^  (11-28) 

p(2)  ^  1(-L.y2)^a)  (II=.29) 

4^^  =  -[(4j^/6)  -  P^]  +  h{>^  +  ih^^^n^^h^  (11=30) 

(p^2)  ^  ^(2)  „  ^(2)^^  ^  „(^U)  ^  ,(!),(!))  •  (ii„3i) 

(2)    (2)  (2)    (2) 

Th«5  second  order  terms  (u   j,  v    p  p    s  h    )  and  the  integration 

variable  (with  respect  to  y)  F(^j,t')  can  be  eliminated  among  (11-27  to 

31)  resulting  in  a  single  equation  for  u    t 


2u^l'  *3u(l»u<l'  *     ^411=   0  (11-32) 


Except  for  the  factor  (l/3)  multiplyiig  the  third  derivative  term  and 
which  can  be  eliminated  by  re  normalizing  the  space  and  time  variables^, 
(11-32)  is  the  same  equation  as  (1^7)    which  describes  the  asymptotic 
plasma  motiono   Aside  from  minor  differences,  (II-32)  was  derived  many 
years  ago  by  Korteweg  and  de  Vries"(l895)  following  an  approximation 
procedure  of  Rayleigh^  our  procedure  is  more  direct  and  indicates  more 
clearly  the  asymptotic  nature  of  ( 11-32  )o   Like  ( Li-7 )  for  plasma  waves, 

the  third  order  non-linear  partial  differential  equation  (11-32) 

1/2 
describes  the  motion  of  water  waves  traveling  with  velocity  (gh  )  ' 


Keller's  (19J4-8)  results  are  also  a  rederivation  of  the  small -amplitude 
equilibrium  solutions  by  scale  transformation  with  respect  to  a  small 
parameter o 
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in  the  positive  x  direction  near  and  including  the  equilibrium  state,, 
A  related  fotorth  order  differential  equation  which  approximately 
describes  waves  traveling  in  both  positive  and  negative  x  directions 
due  to  Boussinesq  has  been  rederived  recently  by  Ursell  (1953)« 
(11-32)  can  be  derived  directly  from  this  fourth  order  equation  by 
applying  the  scale  transformations  (11-8,9)  and  letting  e  — >  Oo 
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